Generation of a wave packet tailored to efficient free space excitation of a single atom 
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We demonstrate the generation of an optical dipole wave suitable for the process of efficiently 
coupling single quanta of light and matter in free space. We employ a parabolic mirror for the 
conversion of a transverse beam mode to a focused dipole wave and show the required spatial and 
temporal shaping of the mode incident onto the mirror. The results include a proof of principle 
correction of the parabolic mirror's aberrations. For the application of exciting an atom with a 
single photon pulse we demonstrate the creation of a suitable temporal pulse envelope. We infer 
coupling strengths of 89% and success probabilities of up to 87% for the application of exciting a 
single atom for the current experimental parameters. 
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I. INTRODUCTION 

In the last decades light-matter interaction has been 
studied in a large variety of experiments not only in many 
body systems such as atomic ensembles, where coherent 
laser pulses are transferred and stored to a collective ex- 
citation [1, 2 , but also on the level of single matter parti- 
cles. Efficient coupling of single atoms to single photons 
has been demonstrated by locating an atom inside an 
optical cavity [3]-[6]. Such efficient coupling is manda- 
tory in many proposals for quantum information proto- 
cols which rely on quantum memories and processors, i.e., 
on the efficient interaction of a quantum state of light 
with the internal quantum states of matter. Light is a 
fast and robust way to transmit quantum information, 
but the transfer over long distances requires quantum 
repeaters based on such quantum memories ^TnS^. The 
rate of information transfer through the quantum chan- 
nel depends on the process efficiency of both storage and 
read-out. Such efficient coupling can be achieved using 
a cavity based scheme and accepting the corresponding 
modification of the mode spectrum from continuous to 
discrete. 

An alternative - and quite natural - approach is the 
coupling in free space. In this approach, the direct cou- 
pling of light to matter has to be maximized instead of 
the coupling to a cavity. Besides their potential use in 
various applications, efficient free space coupling schemes 
offer large design freedom. They enable the investigation 
of light-matter interaction without modifying the den- 
sity of states of the electromagnetic field continuum. Al- 
ternatively one can study the effect of, e.g., a modified 
quantum vacuum [10) . 
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Recently, one could witness increased efforts devoted 
to the topic of efficient free space coupling. When dis- 
cussing the corresponding advances one can clearly dis- 
tinguish between two different processes: One is the elas- 
tic scattering of a laser beam by an interacting particle 
for which no substantial excitation of the matter system 
occurs, i.e., the probability to find the matter system in 
its excited state is low throughout the interaction pro- 
cess. The efficient elastic scattering manifests itself in 
the attenuation of a weak coherent laser beam resonant 
with a transition of the interacting particle, as has been 
shown for single molecules [11 , single quantum dots [12], 
single atoms [13] or single ions [M] . Another hallmark of 
efficient elastic scattering is the phase shift of a weak laser 
beam detuned from resonance, which was demonstrated 
recently for the free space scenario [HI [16] . Correspond- 
ing theoretical contributions can be found in Refs. [ 17i - 
21 . 

The other process treated in recent literature is the 
process of efficient absorption of a photon by a single 
atom, i.e., the matter system is to be brought to its ex- 
cited state with high probability. Note that this should 
not be confused with a 7r-pulse. While a 7r-pulse - typ- 
ically containing many photons - transfers the atom ef- 
ficiently to the excited state the pulse cannot be fully 
absorbed by a single atom. Furthermore, the rotation 
angle of the Bloch vector (e.g. tt) is given by time in- 
tegrating the electric field envelope, while a single pho- 
ton Fock state has a zero electric field expectation value. 
Theoretical and conceptual contributions devoted to the 
absorption of single photons can be found in Refs. [22]- 
27 . A recent experimental demonstration of the absorp- 
tion of a single photon by a single atom was given by 
Piro et al [28 , whereas an early demonstration of re- 
duced transmission of a bright laser beam by a single ion 
was performed by Wineland et al [29 . 

It is common to the two processes of elastic scattering 
and absorption of a photon that their efficiencies are pro- 
portional to the extent to which the incident light mode 
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resembles an electric dipole wave [T9l [2Q l [22 l l23 l [25]. 

More precisely, the spatial properties of the incident light 
mode should match the ones of the wave that would 
emerge from the dipole moment of the interacting parti- 
cle. Here, we show how to experimentally generate such 
a dipole wave envisioning close to optimum coupling ef- 
ficiencies. As outlined below, we do this by preparing 
a plane light wave that, after refiection off a parabolic 
mirror, is transformed into a dipole wave [23l[24]. An es- 
sential property of this setup is that the parabolic mirror 
covers nearly the complete solid angle surrounding the 
dipolar absorber. 

For the specific aim of perfectly exciting a single two- 
level atom with a single photon one may take guidance 
from the spontaneous emission process: The incident 
photon must resemble the time reversed version of a 
spontaneously emitted photon [22H24] . This includes the 
dipole properties discussed above. Furthermore, time re- 
versal arguments demand that the temporal envelope of 
the incident photon is exponentially increasing with a 
time constant equal to the lifetime of the atom's excited 
state [24 . This heuristic approach has been confirmed 
by detailed calculations [26 . 

There exist several methods for temporal shaping of 
single photons in general. One approach is to shape 
the photon during its creation process, for instance by 
temporal modulation of the intensity of the control laser 
in schemes based on electromagnetically induced trans- 
parency [30 or by modulation of the pump beams driv- 
ing the process of four wave mixing ^31j or parametric 
down- conversion (PDC) |32]. Also the shaping of single 
photons by spectrally filtering one of two PDC photons 
has been demonstrated [33 . In the case of single pho- 
ton generation with a single emitter located inside an 
optical resonator [ST-'ST the temporal modulation of the 
pump laser intensity changes the emitted single photon 
waveform. Furthermore a heralded single photon can 
also be modulated directly using electro-optic modula- 
tion HEl]. 

However, with the mentioned techniques it is challeng- 
ing to generate a single photon Fock state with an expo- 
nentially increasing temporal envelop, especially at wave- 
lengths in the ultraviolet spectral range (which are of in- 
terest in this paper, see below). Instead we use coherent 
laser pulses which are highly attenuated to less than one 
photon per pulse. This procedure is often used as a sub- 
stitute for generating single photon Fock states [40^. For 
the sake of simplicity, we will use the term 'single photon' 
also for strongly attenuated coherent states, keeping in 
mind that all excitation probabilities mentioned in this 
paper are strictly valid only for true single photon wave 
packets. Nevertheless, if the amplitude of the coherent 
state is low enough, contributions from Fock states with 
photon numbers larger than one are negligible. It is then 
a good approximation to apply the excitation probabil- 
ities for true single photon states to that instances of 
the coherent state wave packets that would be projected 
onto a single photon state upon measuring the photon 
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TABLE 1. Overview of the dipole transitions treated ex- 
emplarily in this paper. Lifetime data are taken from 
Refs. gl] [42]. The labels Tl and T2 will be used in the 
text for the sake of simplicity. 



number. 

In what follows, we present and characterize experi- 
mental techniques to produce the optimal optical mode 
for perfectly coupling light with an atomic dipole, espe- 
cially for the case of exciting a single atom with a single 
photon. The paper is organized as follows: In the next 
section we motivate a measure for the coupling strength 
between a single atom (or any other matter system) and 
light and relate this measure to the absorption efficiency. 



In Sec. HI the experimental generation of a radially po- 
larized Laguerre-Gauss mode is shown. We present the 
characterization of this beam and calculate the overlap 
of the experimental data with the optimal dipole wave. 
Another property of dipole waves is their uniform phase. 
Therefore, it is essential to compensate for phase aber- 
rations of the device that focuses the incident light to- 
wards the matter particle. This topic is treated in Sec. 
IV} Section [V| is devoted to the shaping of the tempo- 
ral properties of a weak coherent state with the aim of 
maximizing absorption. This aspect as well as the spatial 
mode shaping is exemplified for two specific linear dipole 
transitions of singly and doubly ionized Ytterbium, but 
the corresponding schemes can be adapted to any dipole 
transition. The properties of our exemplary transitions 
are listed in table [T[ Finally, in the last section we esti- 
mate the coupling efficiencies and absorption probabili- 
ties that would be achieved when interfacing the gener- 
ated light waves with the corresponding single ions. 



II. DEFINITION OF THE FREE SPACE 
COUPLING STRENGTH 

The coupling constant of the interaction of light and 
matter in the dipole approximation is given by the Rabi 
frequency g ~ jl ■ E/fi in free space as well as in an 
optical resonator, where jl = ji ■ is the atomic dipole 
moment and E the electric field. Depending on whether 
the treatment is semi-classical or fully quantized, jl alone 
or both jl and E are operators, respectively. In either case 
a large coupling constant requires a large portion of the 
incident electric field to be parallel to the atomic dipole. 

In cavity quantum electro-dynamics the coupling 
strength is proportional to . Here, we adopt this mea- 
sure (cf. Ref. [24 ) and define the free space coupling 
strength G as the electric energy density in the focus 
that corresponds to the field component parallel to the 
atomic dipole and normalize it to the maximum possible 
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one for a given input power: 



(1) 



In the focus, the electric field component along a spe- 
cific dipole moment is proportional to the square root 
of the solid angle covered by the mediating focusing op- 
tics Q.^ and to the overlap 77 of the incident field with 
the dipole field distribution [43 : E - \AV ' V- 

It should be noted that Vt^ is the solid angle weighted 
with the angular dipole radiation pattern D{'d): Vt^ = 
/ D{'d) sini^di^dcp. The overlap 
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is calculated integrating over the solid angle covered by 
the focusing optics [43]. The maximum values of these 
quantities are ^1"^^^ = 87r/3 and rf'^'^ = 1. Therefore, 
the free space coupling strength is written as 
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It reaches unity only if the full solid angle is covered by 
the focusing optics and the mode matching is perfect. 
Likewise the maximum coupling strength for half solid 
angle illumination is G = 0.5. These values correspond 
to the 'scattering ratios' of four and two, respectively, as 
used by Tey et al [21 and Zumofen et al [20 . 

G can be interpreted as a geometric quantity that is 
determined by the focusing setup and the spatial prop- 
erties of the incident light mode. For a classical field or 
a coherent state wave packet a/G scales the field ampli- 
tude at the focus (or its expectation value, respectively). 
If the incident light field is a Fock state, a/G scales the 
magnitude of the uncertainty of the incident field at the 
focus. 

While G = 1 is necessary for reaching the maximum 
effect in all sorts of free space interaction (phase shifts, 
scattered power, absorption, etc.) it is not sufficient for 
reaching unit probability during the absorption of a sin- 
gle photon by a single atom. As mentioned above, the 
temporal profile of the single photon wave packet plays 
a decisive role. The optimum effect, i.e., 100% exci- 
tation of a two-level atom can only be obtained if the 
temporal profile matches an increasing exponential with 
proper time constant [24, 26, 27 . To quantify deviations 
from this profile we introduce the temporal overlap rjt of 
the incident field amplitude Einc{t) with the amplitude 
E[dea\{t) = • 0{—t) of the ideal pulse shape: 
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0{t) is the Heaviside function and F the spontaneous 
emission rate of the two-level system to be excited. This 



overlap can be considered as a generalization of Eq. (8) 
of Ref. [26 which treats the case of an exponentially in- 
creasing pulse amplitude. It is assumed that the incident 
signal is shifted in time by an amount that maximizes 
rjt. Since the respective equation in Ref. [26j determines 
the probability amplitude to find an atom in the excited 
state, we take 77^ as a multiplicative factor determining 
the absorption probability 
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The prospective coupling efficiencies G and absorp- 
tions probabilities Pa based on the experimentally 



achieved overlaps will be given in Sec. VI 



III. 



GENERATION OF THE OPTIMUM 
SPATIAL MODE 



A dipole wave incident from nearly the full solid angle 
is generated by reffection of a collimated laser beam from 
a deep parabolic mirror. The ideal beam shape incident 
on the parabolic mirror is the time reversed version of a 
dipole wave that emerges from the focus of the mirror and 
is reffected at the parabolic surface. After reffection off 
the parabolic mirror the ffeld vector Ej^ of light emitted 
by a linear dipole located at the focal point and with 
the quantization axis along the mirror axis is described 
by IM1I13] 
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1 
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where / is the focal length of the mirror, er the unit vec- 
tor in radial direction, ^0 is an amplitude factor and a 
factor exp(z^r) has been omitted. This ffeld distribution 
has to be generated by an optical element from an in- 
coming plane wave. The amplitude distribution of this 
ideal mode can in principle be created using diffractive 
optical elements. However, there are several drawbacks 
to this method p3 . Instead we generate a radially po- 
larized doughnut beam which can have a large overlap 
with the ideal mode and can be created without using 
diffractive optics. 

The ffeld vector ^rpd of a radially polarized doughnut 
mode can be expressed as 



rpd 
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The overlap of the radially polarized doughnut mode with 
the ideal distribution is optimized by adjusting the beam 
radius w with the maximum achievable overlap depend- 
ing on the magnitude of the weighted solid angle (i.e., 
the depth of the mirror with respect to its focal length or 
the half opening angle) [43 . Our parabolic mirrors have a 



4 



Laser 





Temporal shaping 
1 trigger h AWG -i 


1 — 




1 1 « 

AOIVI 



SIVIF 



Spatial shaping 




Aberration correction 


segmented P'^^oie V 
HWP 




i 13 




phase plate PM 



FIG. 1. Experimental setup. The temporal profile is formed 
using an arbitrary waveform generator (AWG) which deter- 
mines the temporal envelope of the light power scattered into 
the first diffraction order of an acousto optical modulator 
(AOM). The output of the AOM is filtered by a single mode 
fiber (SMF). Linearly polarized light is sent through a seg- 
mented half wave plate (HWP). Subsequently the produced 
doughnut mode is filtered by a pinhole aperture and expanded 
to appropriate dimensions. The aberrations of the parabolic 
mirror (PM) are compensated using a phase plate. 



focal length of / =2.1 mm and a front aperture radius of 
10 mm. This corresponds to a length of 5.67 ■ f (aperture 
radius Tmax/ f = 4.76, half opening angle of 134°) and 
a weighted solid angle for a linear dipole of 1^^=0.94—. 
The maximum possible overlap for this geometric param- 
eters is 77=0.982 for a beam radius w = 2.26 • / [43^. This 
limits the experimentally achievable coupling strength G 
and absorption probability to 0.91. A bore of 0.75 mm 
radius concentric to the optical axis of the mirror (cf. Sec. 
IV) induces only negligible changes to and 77. 

There are different methods for producing a radially 
polarized doughnut beam. Examples are the use of sub- 
wavelength dielectric gratings ^44l [45], nematic liquid 
crystals [46] and the coherent addition of the two modes 
TEMoi and TEMio inside laser cavities [47| or interfer- 
ometers [Ml ^3- We use a segmented half- wave plate 
similar to the one described in Refs. [50l [51] to gener- 
ate the radially polarized doughnut beam. It consists 
of 71 = 8 segments of a low order half-wave plate suit- 
able for both wavelengths (370 nm and 252 nm). The 
experimental setup to generate the doughnut mode is 
shown in Fig. [ij A continuous wave laser beam at the 
wavelength of interest first passes the setup for tempo- 
ral pulse shaping (cf. Sec. |V|). It is then sent through 
a polarization maintaining single mode fiber (370 nm 
beam) or a filtering telescope with a pinhole aperture 
(252 nm beam), respectively. The resulting linearly po- 
larized Gaussian beam is centered onto the segmented 
half-wave plate where its polarization is rotated in every 
segment by an angle of 27r/7i, obtaining a full 27r rota- 
tion of the electrical field vector over the whole beam 
cross section. After propagation through the segmented 



half-wave plate the beam contains the radially polarized 
doughnut mode and modes of higher order. The latter 
ones are filtered out using a pinhole aperture in the fo- 
cal plane of a Keplerian telescope. With this setup the 
doughnut mode is generated with a power efficiency of 
70%. 

The polarization properties of the generated doughnut 
beam are determined performing spatially resolved mea- 
surements of the Stokes parameters based on the method 
described in Ref. [5T . The profile of a beam sent through 
a rotatable quarter-wave plate and a non rotatable lin- 
ear polarizer is measured with a CMOS camera. Due 
to the technical limitations of the CMOS camera the 
Stokes measurements are not performed with single pho- 
ton pulses but with continuous wave laser light. 

Once the Stokes parameters are measured the total 
intensity 5*0 and the characteristic angles of the polariza- 
tion ellipse are calculated for every camera pixel. Since 
the CMOS camera chip used in the Stokes measurements 
is much smaller than the front aperture of the parabolic 
mirror, the measurements have been performed before 
the final magnification step. Nevertheless, the data are 
presented in a coordinate system scaled such that the 
beam radius of the measured doughnut mode matches 
the optimum beam radius predicted by theory. The re- 
sults are shown in Figs. |2] and [3] for 370 nm and 252 nm, 
respectively. 

The orientation angle G [0, tt] describes the rota- 
tion of the main axis of the polarization ellipse with re- 
spect to the X-axis. As shown in Figs. ^) and|3]3), at 
both wavelengths covers the full range of to tt in 
the lower and the upper half plane of the beam cross 
section. The value of 7/^ changes continuously with rota- 
tion around the beam center. The relative phase of the 
electric field is the same in both half planes, otherwise 
a defect line would be observed in at the y = line 
where ip jumps from zero to tt. Hence, the orientation of 
the polarization vector is radial on the majority of the 
beam cross section. Deviations from this behavior occur 
mainly for radii ^ 3.5 • /, i.e., in regions of low inten- 
sity. We attribute these deviations to measurement noise 
originating from the discretization of the camera pixels' 
intensity values. Furthermore, measurement noise is par- 
ticularly cumbersome in the vicinity of the x-axis, since 
there it leads to fluctuations between values of i/j = 
and ip = 7T. When calculating the overlap with an ideal, 
noise-free state of polarization (performed below) these 
fluctuations lead to pixels on which the scalar product 
of the field vectors becomes negative. For pixels with a 
large value of 5*0 this results in a measurement-induced 
reduction of the overlap integral. 

The second angle describing the polarization ellipse is 
the ellipticity angle x ^ [~'7r/4, 7r/4] giving the amount 
of elliptical polarization. Figs. ^) and [sj^) show the 
ellipticity angles over the cross section of the generated 
doughnut beams. On the majority of the beam cross 
section the ellipticity angle is rather low, i.e., the local 
state of polarization can be regarded as linear with a 
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FIG. 2. Characterization of the radially polarized doughnut beam at 370 nm (transition Tl): (a) Stokes parameter So (total 
intensity), (b) Orientation angle ip. (c) Ellipticity angle x- The shown data points cover the aperture of the parabolic mirror. 
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FIG. 3. Characterization of the radially polarized doughnut beam at 252 nm (transition T2): (a) Stokes parameter 5*0 (total 
intensity), (b) Orientation angle ip. (c) Ellipticity angle x- The shown data points cover the aperture of the parabolic mirror. 



radial orientation. 

Next we calculate the overlap of the ideal electric 
dipole field and the generated radially polarized dough- 
nut mode in the entrance plane of the focusing parabolic 
mirror, assuming perfect reflection at the mirror surface. 
This gives the same value as calculating the mode over- 
lap on a sphere close to the focal point, as described in 
Ref. [32. The integration is performed over the mir- 
ror aperture, i.e., only pixels lying within the boundary 
r < Tmax = 4.76 • / are contributing. The spatial field 
overlap rj of the measured radially polarized doughnut 
beam with the ideal dipole distribution reflected at the 
mirror is calculated by taking the normalized scalar prod- 
uct of the electric fields 

^J S\Eme.s?rdrd^J\E^\^rdrdc^ 



where Cmeas is the normalized measured polarization vec- 
tor constructed from the angles ijj and x- ^meas is calcu- 
lated by taking the square root of the measured intensity 
distribution (the Stokes parameter 5'o) and is given 
by Eq. ([6|. Equation (|8| is the equivalent to Eq. 

For the mode driving transition Tl we calculate the 
field overlap of the generated radially polarized doughnut 
mode with the ideal dipole wave to be 0.963. If we assume 
perfect orientation angles for the measured beam, i.e., we 
disregard the discretization effects discussed above, we 
find 77370 = 0.979. The overlap of the generated dough- 
nut mode for transition T2 with the ideal dipole mode is 
calculated to be 0.939. Assuming ideal orientation angles 
we arrive at 77252 = 0.975. 

Our parabolic mirrors are made from aluminum. One 
might wonder whether the radially varying amplitude of 
the reflection coefficient, which is due to the radially vary- 
ing angle of incidence, might influence the obtainable 
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overlaps (see also Sec. [TV] ). Accounting for the radial 
variation of the reflectivity in calculating the overlap of 
the modes given by Eqs. (|6| and ^ leads to a slight 
increase by 0.03% in comparison to the maximum ob- 
tainable mode overlap for constant reflectivity. 

Therefore, we conclude that the obtainable overlaps 
are not influenced by the material properties of aluminum 
in a significant way. 



IV. ABERRATION CORRECTION 

An essential prerequisite for generating dipole modes is 
an aberration free phase front: The phase of the focused 
wave should be uniform on a spherical surface enclosing 
the focal point. Such an ideal spherical wave can only 
be generated if the parabolic mirror is free of surface de- 
viations. Since any realistic device has aberrations (see, 
e.g., Refs. [531455] ). these need to be characterized and 
corrected for. Ideally the wavefront incident onto the 
parabolic mirror can be shaped such as to create an aber- 
ration free inward moving spherical wave after reflection 
off the mirror. 

One method of measuring the imperfections of a 
parabolic mirror is to scan its surface using a profilome- 
ter |55] . But the huge relative depth of our parabola 
makes it hardly accessible for profilometer heads. We 
measure the shape deviations by interferometry with 
nanometer accuracy [53 . Let us briefly recall the core 
features of the interferometer setup, details can be found 
in Ref. [53] . We use an interferometer of the Fizeau type 
at the standard helium-neon laser wavelength XueNe = 
632.8 nm. A reference wave is produced at a wedged 
fused silica plate with surface flatness < A/100. A small 
steel (or glass) sphere located concentric to the focus of 
the parabolic mirror is acting as a null element. Any light 
ray traveling parallel to the optical axis of the parabola is 
reflected towards the sphere, reflected back onto its path 
of incidence and interferometrically overlapped with the 
reference beam. Owing to this double pass geometry, we 
measure twice the phase aberration introduced by the 
mirror. Possible aberrations of the sphere are eliminated 
by averaging over different spheres and different rota- 
tional orientations of each sphere relative to the parabolic 
mirror. The deviations from a perfect parabola obtained 
after removal of misalignment aberrations are finally rep- 
resented by Zernike polynomials. 

Our parabolic mirrors are made of diamond turned alu- 
minum and are produced at the Fraunhofer Institute for 
Applied Optics and Precision Engineering at Jena (see 
Sec. Ill for the geometric parameters). A central on axis 



bore of 1.5 mm diameter is used for locating the sphere 
or any other object such as an ion [56 in the focal point. 

The interferometric measurements are performed with 
a radially polarized doughnut beam centered on the op- 
tical axis of the parabolic mirror. In this particular case 
light is polarized parallel to the plane of incidence at ev- 
ery spot on the surface of the parabolic reflector. 



Figure [4pi) shows the result of such a measurement for 
one of our parabolic mirrors. The peak-to-valley (PV) 
value of the aberrations is 0.53 waves at AneNe- The root 
mean square (RMS) value is 0.09 waves at AneNe- The 
data shown correspond to maximum deviations of about 
130 nm from the ideal parabolic shape [53]. This is of the 
same order of magnitude as reported for the parabolic 
mirror used in Ref. [54 . 

Based on these interferometry results we have pro- 
duced a phase plate which nominally imprints the phase 
conjugate of the measured aberrations onto an otherwise 
flat phase front. The fused silica phase plate is produced 
in house by laser lithography and (selective) reactive ion 
etching. In order to check the quality of aberration com- 
pensation the phase plate is inserted into the interfero- 
metric setup between the parabolic mirror and the flat 
reference plate. While constantly repeating the inter- 
ferometric measurement, position and orientation of the 
phase plate are optimized. As a first optimization crite- 
rion we take the reduction of the PV value of the aber- 
rations while for the fine adjustment the RMS value is 
a better measure, since this is the quantity determining 
the Strehl ratio [57]. The final outcome of such a mea- 
surement series is depicted in Fig. [4]3). The PV value of 
the aberrations is reduced to 0.45 waves at AneNe- The 
RMS value shrinks to 0.033 waves at AneNe- It should be 
noted that the seemingly small reduction of the PV value 
is due to effects occurring at the boundaries of the phase 
plate and the parabolic mirror. Neglecting the outer 0.5 
mm of the aperture the PV value shrinks to 0.16 waves 
at AHeNe ^ud the RMS value is reduced to 0.023 waves. 

To judge the quality of these improvements we have 
performed simulations of the focal intensities based on a 
generalization of the method by Richards and Wolf [58] . 
We assume a wave front according to the aberrations 
shown in Fig. [4] to be incident on a perfectly reflecting 
parabolic mirror. The incident light wave is a radially po- 
larized doughnut mode (see Eq. ^) with the optimum 
beam radius. The resulting focal intensities are normal- 
ized to the focal intensity obtained for an aberration free 
wave front. This yields the Strehl ratios. At the testing 
wavelength of our interferometer (633 nm), the Strehl 
ratio is improved from 68% to 99% by using the correct- 
ing phase plate. In other words, the combined system 
'parabolic-mirror-plus-phase-plate' is practically capable 
of diffraction limited focusing. 

At the wavelengths for experiments with Ytterbium 
ions the Strehl ratios of the uncorrected mirror are 32% 
for 370 nm and 9% for 252 nm. Close to diffraction lim- 
ited performance seems feasible in these cases if the phase 
plate is designed for the respective wavelength consider- 
ing the corresponding refractive index of fused silica. 

Besides that, a phase plate that perfectly compensates 
the aberrations at the testing wavelength AneNe is also 
rather suitable at the wavelengths of the transitions Tl 
and T2. Assuming a phase plate working perfectly at 
the testing wavelength, we simulated the performance 
of the compound system 'phase-plate-plus-paraboloid' at 
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FIG. 4. (a) Measured phase aberrations in single pass reflection. The resulting aberrations are shown averaged over two 
different spheres and six rotational settings per sphere. The PV value is 0.53 waves, (b) Measured aberrations with a 
compensating phase plate (single experimental run). The PV value is 0.45 waves. In both panels the data result from a Zernike 
fit of degree 10. Adjacent contour lines are separated by 0.05 waves at AneNe- 



the other two wavelengths. Since the relative refractive 
index changes of fused silica at the wavelengths of interest 
are on the order of 1% - 3%, our simulations reveal Strehl 
ratios ranging from 97% to more than 99%, depending on 
the wavelength and the specific parabolic mirror. 

So far we have assumed a phase shift upon reflection 
that does not depend on the distance to the mirror's op- 
tical axis. But due to the complex refractive index of 
aluminum this assumption is generally not justified. The 
radially varying phase shift upon refiection imposes a de- 
focus. Nevertheless we argue that the influence of the 
refractive index of aluminum can be neglected. We calcu- 
late the phase shift upon reflection by taking the respec- 
tive refractive indices from Ref. [59 and determining the 
argument of the complex reflection coefficient for electric 
fields parallel to the plane of incidence. For this we use 
the relations between angle of incidence and radial dis- 
tance to the optical axis given in Ref. [23 . The obtained 
phase shifts are then used in simulations with an oth- 
erwise aberration free paraboloid. For the wavelengths 
considered here the simulations show an axial shift of 
the intensity maximum by less than a tenth of a wave- 
length. The maximum intensity is reduced by amounts 
on the order of 0.2% in the case of transition T2 and even 
less for the other wavelengths. The intensity in the nom- 
inal focus is reduced by 0.5% (1.5%, 2.8%) at the testing 
(transition Tl, transition T2) wavelength. 

The phase shifts induced by the above effects can in 
principle be measured by our interferometric setup. How- 
ever, a substantial part of the treatment of the interfero- 
metric data is the removal of misalignment aberrations of 
the null sphere from the exact focal position [53 . Since 
the monotonous radial variation of the phase shift upon 
refiection is equivalent to a slight displacement of the null 
sphere, this effect is removed in the measurement proce- 
dure. 

We finish this section by noting that the orders of mag- 



nitude of the experimental results shown are generally re- 
produced for other parabolic mirrors, with the shape of 
the aberrations of the uncorrected mirrors differing from 
device to device. 



V. GENERATION OF THE OPTIMUM 
TEMPORAL MODE 

To design the temporal mode of a single photon wave 
packet that would be absorbed efficiently by a two level 
system, we exploit the time reversal symmetry of the 
compound system atom and free space field modes [24] . 
An excited two level atom spontaneously decays into the 
ground state under emission of a photon with a proba- 
bility amplitude that is exponentially decreasing in time 
[60 . If a photon wave packet with such a decreasing enve- 
lope is sent onto a two level system the maximum possible 
excitation probability is calculated to be 54% [26] , 

An excitation probability of 100% is achieved by gener- 
ating a photon wave packet with the time reversed shape 
of the spontaneously emitted photon: an exponentially 
increasing field amplitude with a time constant which 
is twice the lifetime of the excited state. Furthermore 
the ideal pulse length is infinite and has a sharp edge 
at the end. The excitation probability depends not only 
on a good overlap of the time constant and the sharp 
edge at the end of the pulse, but also on the total pulse 
length. In Ref. [26 it is shown theoretically that for pulse 
lengths of five lifetimes 99% excitation probability can be 
achieved. These numbers are valid for single photon Fock 
states. However, as outlined already in Sec. [l| coherent 
state pulses with average photon numbers <Cl are a good 
approximation to single photon wave packets. 

In the following, we discuss the generation of these co- 
herent state wave packets. The experimental apparatus 
for the generation of arbitrarily shaped pulses is shown in 
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Fig. [T] A continuous wave laser beam is cut into pulses 
using an acousto-optic modulator (AOM). The intensity 
diffracted into the first order of the AOM depends on the 
amplitude Uq of the applied radio frequency (RF) signal 
as I {t) ~ sin^([/o(t)). An exponentially increasing wave 
packet is obtained by modulating the amplitude of the 
applied RF signal /7RF(t) = Uo{t) ■ sin(co'RFt) as 



Uo{t) = arcsin ( exp 



t 

27 



(9) 



where cjrf = 27r • 400 MHz is the RF driving frequency 
of the AOM. The modulated RF signal is generated with 
a digital arbitrary waveform generator and amplified be- 
fore being applied to the AOM. The first diffraction or- 
der of the AOM is coupled to a polarization maintaining 
single mode fiber and sent to the setup for the spatial 
beam shaping. The attenuation to a mean photon num- 
ber per pulse of much less than one is done using a set 
of half- wave plates and linear polarizers. The temporal 
profile of these highly attenuated laser pulses is measured 
using a photomultiplier tube with a quantum efficiency 
of 15% and a frequency counter in start-stop configura- 
tion. Since we do not measure a photon in every pulse, 
the start trigger is given by the detection of a photon 
and the stop trigger by the following clock signal acti- 
vating the output of the arbitrary waveform generator. 
We account for this trigger configuration by reversing the 
temporal coordinate of the experimental data. In Fig. |5] 
the measurements of the pulses designed to match the 
lifetimes of transition Tl and T2 are shown. The mean 
photon number per pulse is 0.09. For simplicity reasons 
the 370 nm laser was used for the creation of the pulses 
for both time constants. 

The temporal overlap r]t of the field amplitudes is cal- 
culated according to Eq. Q. The incident field am- 
plitude Ei^c{t) is obtained by taking the square root of 
the photon count statistics and the integration is per- 
formed over the time intervals covered by the experimen- 
tal data shown in Fig. [5] The results of this procedure 
are r]t =0.96 for the pulse with 8 ns time constant and 
r]t =0.99 for the 230 ns pulse. The lower overlap value 
in the case of the 8 ns pulse is attributed to the finite 
build- up/decay time (5 ns) of the optical grating inside 
the AOM crystal. This results in non-negligible devia- 
tions from a sharp pulse end when compared to the time 
constant of the rising exponential. 



VI. DISCUSSION 

With the acquired experimental data at hand we are 
now able to give an estimate of the coupling efficiency 
and absorption probability which would be achieved in 
experiments employing the generated field modes. Ta- 
ble |ll| summarizes the field overlaps of the generated light 
modes in time and in space domain, as well as the Strehl 
ratio. According to Eqs. ([3| and ^ we calculate G 




100 



2000 



FIG. 5. Generated exponentially increasing pulses with a 
mean photon number per pulse of 0.09 photons. The dashed 
line indicates the perfect pulse shapes. The time constant 
(life time) is 8 ns in panel (a) and 230 ns in panel (b). The 
total number of counts are 430 • 10^ and 487 • 10^ respectively. 



Transition 


Tl 


T2 


Spatial overlap r/ 


0.979 


0.975 


Solid angle ^^/^ 


0.94 


0.94 


Strehl ratio 


0.99 


0.99 


Temporal overlap r/t 


0.96 


0.99 


Coupling strength G 


0.892 


0.885 


Absorption probability Pa 


0.812 


0.867 



TABLE II. Overview of the achieved figures of merit. The 
Strehl ratio is taken from the proof of principle experiment at 

AneNe- 



and Pa from these values and let the Strehl ratio en- 
ter as a multiplicative factor accounting for losses caused 
by the mirror imperfections remaining after compensa- 
tion. We thus arrive at prospective coupling strengths of 
G = 0.892 (G = 0.885) for transition Tl (T2). These 
values are only about two percent below the upper limit 
set by the experimental geometry when using doughnut 
modes. The corresponding absorption probabilities are 
= 0.812 (P, = 0.867). 

These absorption probabilities are valid under the as- 
sumption of a closed two- level system (e.g., transition 
T2 of ^^^Yh^^). As an example for a more complex 
level structure we take the transition Tl of ^^^Yb+: 
Both the excited and the ground state are split into two 
Zeeman sub-levels with the magnetic quantum numbers 
m = ±1/2 (we neglect the dipole allowed ^^3/2 
transition due to the small branching ratio of 0.5% [61]). 
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Two linear and two circular dipole transitions are possi- 
ble. The probability for spontaneous decay via a linear 
dipole transition is 1/3 according to the corresponding 
Clebsch-Gordan coefficient and 2/3 via a circular dipole 
transition. Exciting the linear dipole transition leads to 
a maximum possible excitation probability of 1/3, as- 
suming the excitation light mode is ideal. In this sense, 
one has to understand the value Pa = 0.812 given for 
a mode designed for this transition as a relative value, 
i.e., the de facto absorption probability would be Pa/3. 
However, the coupling strength G itself is not affected by 
such considerations. For example, an elastic scattering 
experiment on any of these transitions should be feasi- 
ble with the coupling strength derived for an ideal two- 
level system because no dipole moments are induced for 
transitions other than the driven one (cf. Ref. [62 for a 
calculation applying the level structure discussed above). 

We conclude this paper by emphasizing that the 
method presented here to generate a dipole mode can be 
applied to a dipole transition at any optical wavelength. 



An application to circular dipole transitions is straight- 
forward requiring corresponding amplitude and polariza- 
tion distributions. However, these polarization patterns 
become quite complex - especially when focusing from a 
large solid angle [43]. 

In summary, we have demonstrated the generation of 
modes of the electromagnetic field allowing for efficient 
coupling of single photons and single two-level systems 
in free space using a deep parabolic mirror. The free 
space coupling method described here should compete 
favorably with other coupling schemes. 
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